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Abstract 

We study the median m(x) in the gamma distribution with parameter 
x and show that < m'(x) < 1 for all x > 0. This proves a generalization 
of a conjecture of Chen and Rubin from 1986: The sequence m{n) — n 
decreases for n > 1. We also describe the asymptotic behaviour of m and 
m' at zero and at infinity. 
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1 Introduction 

The gamma distribution with parameter x has density with respect to Lebesgue 
measure on (0, oo) given by e~ l t x ~ x /T(x). We consider the median m{x) of this 
distribution which is defined implicitly as 

or 

(•m(x) i roo 

e'H^ 1 dt = - e-'i*- 1 alt. (1) 



o 2 j 



We of course also have 

'OO -i roc 

e'H x - 1 dt = - e-H x - x alt. (2) 

m{x) 2 Jo 

We show that m is continuous and increasing. This is a consequence of a result 
about general convolution semigroups of probabilities on the positive half-line, 
that we give in Section [51 There we also show that m is real analytic and that 
m satisfies a certain differential equation. 

We shall mainly study m through the function 

x 

(f{x) = log— r-r, x > 0. (3) 
m(x) 
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This function appears if we make the substitution u = \og(x/t) in the relation 
Q. We get: 

/tp(x) _ 1 poo 

e -x(e --Hi) du=~ e~ x( - e du. (4) 
-oo * J — oo 

Chen and Rubin (see jSJ) studied the median of the gamma distribution and 
proved that x — 1/3 < m(x) < x for x > 0. The relation (J3J was also used in [2] 
to establish that <p(x) > 0, or equivalently m(x) < x. (It follows by observing 
that f^ oo e~ x ( e u+u ~) du < J °° e~ x ( e u + u ) du, which is true because sinhn > u 
for u > 0.) 

Chen and Rubin furthermore conjectured that the sequence m(n) — n de- 
creases. This conjecture has recently been verified by Alzer (see (J). He proved 
that m{n + 1) — an decreases for all n > exactly when a > 1 and increases 
exactly when a < m(2) — log 2. Note that m(l) = log 2. 

In this paper we investigate the properties of m as a function on (0, oo). 
Chen and Rubin's conjecture follows from the relation 

m'(x) < 1 for all x > 

and it is our main goal to verify this relation. (See Theorem 13 .31 ) 

We have drawn the graph of m in Figure ^ (The dotted line is given by 
x- 1/3.) 




x 



Figure 1: The graph of m 

In terms of ip Chen and Rubins conjecture takes the form l — xip'(x) < e^ x >. 
As a crucial step towards this result we show that xip(x) decreases and that 

- < xip(x) < log 2 
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for all x > (see Proposition 13 . 5 j) . In terms of m this relation can be rewritten 
as 

xe -io g 2/x < m ^ < xe -i/3x_ 
If we use that e~ a < 1 — a + a 2 /2 for a > then we see that in fact 

. . / 1 1 \ 11 

mix) < x 1 1 =■ = x j 

V ' V 3x 18xV 3 18x 

for all x > 0. This result improves a result of Choi about m(n + 1), see [HI 
Theorem 1], even though it is there claimed to be best possible. 

The asymptotic behaviour of m(x) for x — > and for x — * oo is given in 
Section |1J From this we can deduce that vrS- k \x) > for x close to for each 
k while (— l) k m^ k \x) > for x sufficiently large and k > 2. It is reasonable to 
believe that m"(x) > for all x > 0, i.e. m is convex, but the higher derivatives 
of odd order change sign. 

We also relate our work to papers of Ramanujan, Watson and others on 
Ramanujan's rational approximation to e x , see Section |HJ 

2 Medians of convolution semigroups on the half- 
line 

A family {fJ- x }x>o of probabilities concentrated on [0, oo) is called a convolution 
semigroup if it has the properties 

(i) n x ([0, oo)) = 1 for all x > 0; 

(ii) fi x * fiy = /J, x+y for all x, y > 0; 

(hi) n x — * Sq for x — * in the vague topology (here 6q denotes the Dirac mass 
at zero). 

The gamma distributions {e~ t t x ~ 1 /T(x) dt} x> o is an example of such a convo- 
lution semigroup. 

A probability measure [i on [0, oo) has median m if 

M ([0,m]) = i. 

Of course a probability measure may not have a median and if it exists it may 
not be unique. However, if the measure has density w.r.t. Lebesgue measure on 
[0, oo) then the median exists: this follows from the fact that M — * /i([0, M]) is 
continuous and increases from to 1. If the density is strictly positive almost 
everywhere on [0, oo) then the median is unique. 

Proposition 2.1 Let {[i x }x>o be a convolution semigroup of probabilities on 
[0, oo) having a.e. strictly positive densities w.r.t. Lebesgue measure on [0, oo). 
Then the median m(x) of fi x is a continuous and strictly increasing function 
on (0, oo). 



3 



Proof. We write dfx x (t) = g x (t)dt. As we have seen above, m exists as a 
function on the positive half-line. We now show that m is increasing. Let x > 
and h > 0. Then 



/ g x+ h(t)dt = 
Jo 



f 

Jo 

f 

Jo 

f 

Jo 



< 



9x *gh(t)dt 

/ 9h(t - s)g x (s)dsdt 
Jo 

/m(x) 
g h (t - s)dtg x (s)ds 

9h{t) dt g x (s) ds 

poo 

/ 9h{t) dt g x {s) ds 
Jo 

rm(x) 

g x (s)ds = -. 



m{x)- 



Since m(x + h) by definition satisfies 

«m(x+h) 



f 

Jo 



9x+h{t)dt = -, 



we must therefore have m(x) < m(x + h). This shows that m is strictly increas- 
ing. 

Concerning the continuity we first notice that for any A > 0, 



9x +h(t)dt -> / g xo (t)dt 



as h — > 0. In fact, using the same computation as above we get 

r-A r-A 

\ 9x a +h{t)dt = I g X0 *gh(t)dt 
Jo Jo 

A pA-s 

\ 9h(t) dt g xo (s) ds 
o Jo 

r-A r-A-s 

-> dS (t) g XQ (s) ds 

Jo Jo 

= / 9x (s)ds, 
Jo 

as h — > 0, because of the vague convergence. (All measures have the same total 
mass, so vague convergence implies weak convergence.) Now let xo > and let 
e > be given. For A = to(xq) + e we have, as h — ► 0+, 



/ 9x +h(t) dt-> 

Jo Jo 



m(xo)+e 



g xo {t)dt > -, 
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so that m{xQ + h) < m(xo) + e for all sufficiently small and positive h. Similarly 
we find m(xo — h) > m(xo) — e, and this shows that m is continuous. □ 
We now specialize to consider the medians m(x) of the gamma distributions. 
We notice 

Proposition 2.2 The median m(x) of the gamma distributions is real analytic. 
Proof. We consider the C 1 -function 

F(x,y) = [ V e-H^-^-dt, x,y>0. 
Jo r W 

The median is implicitly defined as F(x, m{x)) = 1/2. The fact that the contin- 
uous function m is C 1 follows from the implicit function theorem, which yields 
the following differential equation for m: 

1 pm(x) 

e- m ^m(x) x - 1 m'(x) = -T'{x)- (logt) e'H^ 1 dt, (5) 

2 Jo 

which shows that m satisfies a differential equation of the form 

m'(x) = G(x,m(x)), 

with G(x,y) being real analytic for x,y > 0. Therefore m is real analytic. □ 

Remark 2.3 From |3p it seems difficult to deduce monotonicity properties of 
m; e.g. it is not at all clear that m'{x) > 0, which we know from Proposition 
\2.1\ For another derivation ofm'(x) > see Provosition VS . 51 

3 Uniform results 

In this section we prove the generalized version of Chen and Rubin's conjecture 
(Theorem ESI). 

Proposition 13.21 below is the key to our results. Before stating it we need 
some notation. 

We consider the function f(x) = e~ x +x. It is easily seen that f(x) decreases 
for x < and then increases for x > 0. Hence / has an inverse on (— oo,0], 
which we call u, and an inverse on [0, oo), which we call v. The function u is 
defined for t > 1 as u(t) < and 

e~ u{t) + u(t) = t, t > 1; 
the function v as v(t) > and 

e -«(*) + v (t) =t, t>l. 
The following function £ plays an important role: 

m = tt^w + rd^y = «'(*) + v 'v- ( 6 ) 

The function £(t) is defined for t > 1 by ©, but for t = 1 this expression yields 
oo — oo and a closer study is necessary. The following result holds and will be 
proved in Section 
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Proposition 3.1 The function ^(t) defined for t > 1 by {j^ has a holomorphic 
extension to the cut plane C \ {x ± 2m \x > 1}. The following holds: 

(i) lim^-oo^t) = 0, lim^oo^i) = 1, £(t) is increasing; 

(ii) lim^oo ^ n \t) = forn> 1; 

(Hi) is concave for t £ [l,oo); 

(iv) £(i) = 2/3, £'(1) = 8/135 and f" (1) = -16/2835. 

('vj ^4Z/ points on {x ± 27ri | x > 1} are singular points for £ and i/ie series 



fc=0 



diverges for all x. 
In Figure [21 we have shown the graph of for i > 1. 




Figure 2: The graph of £ 
Proposition 3.2 We have 

ftp(x) POO 

2 J e -*(e-"+«) = y £(i)e _:!;t dt, (7) 
where the function £ is defined in 
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Proof. From (J3J we get 

r<p(x) roo rO 

2 / e ~ x{e ~ u+u) du= / e~ x{e ~ u+u) du - / e~ x{e ~ u+u) du. 

JO JO J-oo 

In the first integral on the right hand side we make the substitution u = v(t) 
and in the second the substitution u = u(t), where u(t) and v(t) denote the 
functions above. In this way we get 



dt 



/oo 
i{t)e~ xt dt. 

□ 

It is easy to see that in terms of the function (p the relation m'(x) < 1 takes 
the form 

1-VW < e^ x \ 



Theorem 3.3 We have 

1 - x<p'(x) < e^ x) 

for all x > 0. 

Before proving this result, we give the following Lemma. 
Lemma 3.4 The relation 

pxip(x) n poo 

Jo 3 Ji 

holds for all x > 0. 

Proof. In the left hand side of the relation in Proposition 13.21 we make the 
substitution s = xu and get in this way 



rx(p(x) poo 

J e - s - xe ' 8/x ds = J £(t)e- xt dt. 



2 r x ^ _.___./. . r°° 

x Jo 

On the right hand side we perform integration by parts, and this gives 



rxip(x) roo 

2 J e~ s - xe 1 ds = J i{t)xe~ xt dt 



oo 

= e~ x i{l) + J £'{t)e- xt dt. 

Here £(1) = 2/3 (see Proposition 13.11) and the assertion of the lemma now 
follows by multiplication by e x . □ 
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Proposition 3.5 The function x — > x(p(x) decreases for x > anc? we /iai>e 

lim x</?(aj) = log 2, 
i— >o+ 

lim xy(aj) = — . 

x^oo 3 

In particular <p is decreasing from oo to (and m is increasing). (The graph of 
xip{x) is shown in Figure^) 




Figure 3: The graph of x<p(x) 
Proof. From Lemma 13.41 we find by differentiation, that 

e -»(v(*)-l+e-*->) ( X(p ( x ))' = 
rxip(x) 



— / te~ xt C{t + l)dt. 

Since 1 - (1 + a)e~ a > for a > and £'(t -f 1) > for t > (see Proposition 
13.1(1 we have 

= J e -^ e -(i- e - s/a; ) (l - (l + ~) e"'/*) > 

and 

S(x) = - / ie-^'O + l)^ > 0. 
2 Jo 



Therefore we have 

(x<p(x))' = -e x{ -^- 1+e ' v{x) \A{x) + B(x)), (8) 

which is a negative quantity, so xip(x) decreases. 

Furthermore, since e - s e x ( 1 ~ e s/x ) > e _s , we get from Lemma 13.41 

rxip(x) n /-oo o 

2/ e~ s ds<-+ £( t )dt = -+ lim£(x)-£(l) = l, 

since lim^oo = 1. (See Proposition 13.11 ) From this we obtain that xip(x) 
is bounded by log 2 for all x > 0. 

If we let I = lim x ^o + x<p(x) then by the dominated convergence theorem we 
find 



rl ry /"OO 

;J e - s ds = - + J ?(t)dt = l, 



so that I = log 2. 

We let finally L = lim^^oo xip(x). We get in the same way as before 

2 

1* = -, 

so that L = 1/3. □ 
Proof of Theorem Iff. 51 We use the expressions A(x) and B(x) and the 
relation @ from the proof of Proposition 13.51 We get 

l-V(z) = e xMx) - 1+e ~ v{x)) (A{x)+B(x)) + ^{x) + l 
< e M*) (A(x) +B(x)) + tp(x) + 1. 

Now, using that 1 — e~ a < a and 1 — (1 + a)e~ a < a 2 /2 for a > 0, we find 

. , - / s 2 \ x(^(x) 3 



pxip(x) / 2 \ 
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Furthermore, since £'(t + 1) < £'(1) = 8/135 for i > (see again Proposition 

EB, 

This gives 

i - VM < i + *.) + ^" + ^g^) 

It is easily seen that the function 



p{u) = e" - + 



4 \ 



6 ' 135m 2 J 

attains its maximum on [1/3, log 2] at u = 1/3 with value 

r V 3 J 90 2 



9 



We obtain from this the relation 

1 - x<ff{x) < 1 + <p{x) + ^{xf < e^ x) . 

□ 

Remark 3.6 The difference between e^ x ' and 1 + (f(x) + ^cp(x) 2 is 0((p(x) 3 ), 
that is (by Proposition^^) the difference is 0(x~ 3 ). Hence for large x the 
difference is very small, reflecting the fact that m'{x) is very close to 1. For 
x close to 0, the difference is large, reflecting the fact that m'(x) approaches 
rapidly. See Proposition \4- 1\ 

4 Asymptotic results 

Here we describe the behaviour of m and m! at zero and at infinity. We first 
investigate m near zero. We use f(x) ~ g(x) for x — » to denote that 

\im f(x)/g{x) = 1. 
a;— >0 



Proposition 4.1 We have m{x) x — > 1/2 as x — > and 

m(x) ~ e" 7 2 _1 / a ' as x — » 0. 

Fwri/iermore, 

2-l/a. 

m'(x) ~ (log2)e 7 — ^— as x — > 0. 

x z 

Proof. For the function l{x) = log(m(x) :r ) we have l(x) = —x(p(x) + xlogx. 
Since xlogx — > as re — > then — ► — log 2 as x — > 0. Therefore m{x) x — > 
1/2 as s -► 0. 

By definition of m(x) and by the functional equation of the gamma function 
we have 

rm(x) i 
/ e - t xt x - 1 dt = -T(x + l). 
Jo 2 

We perform integration by parts on the integral on the left hand side in this 
relation and we thus get 

rm(x) 1 
m (x) x e- m( - x) + / e~H x dt = ~T(x + 1). 
io 2 

We next differentiate this relation and get after some manipulation 

-r'(x + l)- / (logrte-^rft I . 
2 Jo / 

We now use that m(x) x — > 1/2 as x — > and we get in this way 

lim ( log mix) H m'(x) ) = T'fl) = -7. (9) 

x-*o \ m(x) J 
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This is the same as l'(x) = (log(m(x) x ))' — ► —7 as x — > 0. Using l'Hospital's 
rule we get 

Z(x) + log2 
> -7 

x 

for x — > 0. Therefore m(x) ~ 2~ 1 / a; e _7 as x — > 0. 
By (0 we get 

x 2 

xlogm(x) H —m'(x) — > 

m(x) 

for x — > and, since (as used before) xlogm(x) — ► — log 2, we get 



m (x) ~ — =- 2 

x^ 



as x — > 0. □ 
At infinity we have 

Proposition 4.2 The functions m and m! have asymptotic expansions at in- 
finity. We have in particular 

. 8 368 . o. 

m (x) = 1 ^ 7? + o(x ) as x —> 00 

v y 405x 2 25515x 3 v ; 

and 

18 184 . _ 2 . 

mix) = x 1 1 tt + o(x as x — * 00. 

1 ; 3 405x 25515x 2 v ' 

Remark 4.3 Choi (^61) found the asymptotic expansion of m(n+l) up to order 
o(n~ 3 ). Higher order expansions of m{n + \) were found in ' L 9l. In the appendix 
we have included higher order expansions of m{x) and <p(x) (and higher order 
derivatives of the function £ at 1). Because of the complexity, the computations 
behind the expansions in the appendix were made using "Maple 9", using the 
same method as described in Lemma \4-4\ 

Lemma 4.4 The functions ip and <p' have asymptotic expansions at infinity. 
For the function ip we have in particular 

. , 1 29 37 . _ 3 , 

G3(X) = 1 7T T + 0(X 

^ V ; 3x 810x 2 25515x 3 v ; 

as x — > 00. 

Proof. We have already seen in Proposition 13 . 51 that ip(x) = 0(x~ 1 ). Therefore 

r<p(x) n (_-[\k k rf(x) 

2/ e -x(n+ e du = 2 { L ± X j ( n + e -«_l)fc ( iu + o(x-( n+1 )). 

fc=0 ' ^° 

On the other hand, by partial integration, 

i(t + l)e~ xt dt = J2 + o{x-^). (10) 

Jo k=o x 
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(We remark here that the sum involving the derivatives of £ at 1 describes the 
asymptotic behaviour of Ramanujans function, see Section HO) From Proposi- 
tion E21 w e thus get 



^) = Eiirr-E L ^ ! r W (« + e--i)*Ai + o(x-(^)). (ii) 

k=0 k=l J0 



From this relation it is possible to deduce that (p(x) has an asymptotic expansion 
of the form 

n 

"CO (12) 

k=l X 

as x — > oo and for any n. (See below.) 

Let us first find the coefficients c\, C2 and C3 in the expansion. If n = we 
get, using £(1) = 2/3, 

<p(x) = — + o(x~ 1 ). 

(This we have already found in Proposition 13.51 ) For n = 1 we get, using 
^(1) = 8/135, 

1 4 f^ {x) 
<p(x) = — — h — — — 7T -j- x I (u + e~ u - l)du + o(x~ 2 ). 



3x 135a; 2 

Here we use u + e~ u — 1 = u 2 /2 + o(u 2 ) for u — > and we get in this way 
<P(x) = ir + T^T + x (^7T- + o{x- 3 )) + o(x- 2 ). 



3x 135x 2 \ 6 

Then we use the result for n = to obtain that x(/?(x) 3 = l/(27x 2 ) + o(x -2 ), 
and if we insert this into the relation above we get 

= h 1 ) ^ + o(x z ) = 1 = + o(x ). 

^ ' 3x V 135 1627 x 2 V ' 3x 810x 2 V 7 

We repeat the argument to obtain the following for n = 2: 

14 8 

= 7T + 



3x 135x 2 2835x 3 

ip(x) 2 /-^(a;) 



+x / (it + e" u - 1) du / (u + e~ u - l) 2 du + o(x~ 3 ), 

where we use u + e~ u — l = u 2 /2 — u 3 /6 + o(u 3 ) and (ii+e -u — l) 2 = u 4 /4+o(u 4 ) 
for u — > and we get 

14 8 

iPW = T" + 



3x 135x 2 2835x 3 
ip(x) 3 ^{x) 4 



+ x ( -Y- - -it + 0(x } 

X 2 / (i?(x) 5 . r.\ . _o 

' + o(x 5 ) + o(x 3 ). 



2 V 20 
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In the term xtp(x) 3 /6 we substitute the expansion of <p(x) for n = 1 and in the 
terms x(p(x) /24 and x 2 <^(x) 5 /40 the expansion of (p(x) for n = 0. In this way 
we get 

1 29 / 8 1 29 1 1 1 1 \ 1 , , s 

mix) = I 1- 1 \-OiX~ ) 

^ V ; 3x 810x 2 V 2835 18 810 24 81 40 243/ x 3 V ; 
1 29 37 

= 1 n 5"+0(x d ). 

3x 810x 2 25515x 3 v ; 

These computations also indicate how to show that there is an asymptotic 
expansion (|12|) for every n > 1. One could use an inductive argument based on 
the relation (Jll[): First of all, the sum 

/ / 2x k+1 

k=0 

contains a term of order l/x n+l . Next the integrand (u + e~ u — l) k is approxi- 
mated by its Taylor polynomial of order n + k, 

n+k 

n+k\ 



l=2k 

and the expansion ip(x) = Y2k=i c k/ x ~ k + o(x~ n ) is then used in the upper limit 
of the integrals in the sum 

Using these approximations it is possible to obtain 

k=0 k=l l=2k \j=l J 

This is an expansion of tp{x) of order n + 1. 

To see that also (p'(x) has an asymptotic expansion we differentiate and 

get 

2^(x)e- x ^ (a;)+e ^ - 2 / e -*(«+« + e -«) du = - ti{t)e~ xt dt. 



Adding Q to the relation above we get after multiplication by e x and a change 
of variable s = t — 1, 

2^(x)e- :E ^( :c ) +e "" W - 1 ) 

= 2 / e -*(«+e- u -l) ( n + e "« _l) du _ S £(l + s ) e -*s d S. 

Jo Jo 

In the first integral we consider again Taylor approximation and in the second 
integral we perform integration by parts. Using also the asymptotic expansion of 
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e x(ip(x)+e <e( x —i) we are g na vjy a \y\ e to see th a t there is an asymptotic expansion 
of <p'(x). The coefficients in this expansion can be identified by integrating the 
expansion and using the known expansion of tp(x), cf 2 [Appendix C]. □ 
Proof of Proposition\4~^ Since (p(x) has an asymptotic expansion, 



m(x) = xe ^ x > = x 



<p(x) -^ST 



l)V(x) 



^ k\ 

k=0 



also has an asymptotic expansion. We have in particular 

m( x ) = Xe -V(*) = x (l- (p( x ) + - + o(x- 

\ 2 6 

We insert in this relation the asymptotic expansion of ip(x) from the lemma 
above. We get, after some computation, 

/ 1 8 184 

mix) = x\l 1 k -\ o + o(x ' 

y ' \ 3x 405x 2 25515x 3 v 

18 184 _ 2 

" X ~ 3 + 40fe + 25515a; 2 + ° {X '■ 

Since (p'(x) has an asymptotic expansion, the same is true for m'(x), since 
m'(x) = (1 — x(p' (x))e~^ x \ The expansion of m(x) can be found by integrating 
the expansion of m'(x), and this gives the desired expansion of m'(x). □. 

Remark 4.5 It is clear that the methods above can be continued so the asymp- 
totic behaviour ofm^(x), k > 2 for x — > and x — > oo can be determined by 
differentiation of the asymptotic formulas for m'(x). 



5 Properties of the auxiliary function £ 

In this section we derive the properties of the function £, stated in Proposition 
13.11 Most of these properties we found surprisingly difficult to establish. Much 
of the difficulty lies in the fact that £ in given as a sum of two terms, where it 
is necessary to control the cancellation between these two terms. 

Throughout this section u and v denote the functions defined in (fB|). We 
begin by considering them separately. The investigation is based on a simple 
lemma of independent interest. We recall that a C°°-function v{t) defined on 
the positive half-line is called completely monotonic if 

(_l)P l/ (p)( t ) > o, for all integers p > 0. 

It is called a Bernstein function if 

v(t) > and (-l) p u^\t) < 0, for all integers p > 1. 

The last conditions can also be expressed that v'{t) is completely monotonic. 
For details about these classes of functions see e.g. 
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Lemma 5.1 Let F be completely monotonic, letc(t) be a positive C°° -function 
for t > and assume that c'(t) = F(c(t)). Then we have: 

For each n > 1 there exists a completely monotonic function F n such that 

(-l) n -\ {n \t) = F n (,(t)). 

In particular c is a Bernstein function. 

Proof. For n = 1 we can use F\ = F. If (— l) n-1 ^( n ) (£) = F n (q(t)) for some 
completely monotonic function F n then 

= -F'Mt)y(t) = -F' n {,(t))F{,(t)). 

Here F n+ \ = —F^F is completely monotonic as a product of two completely 
monotonic fuctions. 

Furthermore we have 

(-l)\ (n \t) = -F n (,(t))<0, 

so q is a Bernstein function. □ 

Proposition 5.2 The functions —u(t + l) and v(t + l) are Bernstein functions 
and 

]imv(t)/t = l, limu(£)/£ = 0. 

t—>oo t— >oo 

Proof. We have v'(t) = F(v(t)) where 



oo 

kx 



e 

k=0 



is completely monotonic. Since v(t) —* oo as t — > oo we have e U W — > and 
hence v(t)/t = 1 - e~ v W /£ -> 1 as t oo. 

For the function it; = — u we have w/(t) = G(w(t)) where 



1 oo 
k=l 



e kx 



is completely monotonic. We have t = e~ u ^ +u(t) > l+u(£) 2 /2 (since u(t) < 0) 
so that —u(t) < \Jl{t — 1) and thus u(t)/t — > as £ — > oo. □ 
Bernstein functions admit integral representations (see e.g. p. 64]) and 
thus we have 

poo 

v(t + l) = t+ / (1 - e~ xt ) dX(x) 
Jo 



and 



/•oo 

-u(£ + l) = / (I - e' xt )da(x), 
Jo 



for some positive measures A and a on (0, oo). Since £+1 — v(t+l) = e 
for £ — ► oo we conclude that A((0, oo)) = 1, hence 



/•oo 

v(t + l) = t + l- / e~ xt d\(x). 
Jo 
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The measure a has infinite total mass because u(t) — > — oo for t — ► oo. 

From these representations we see that = u'(t) + v'(t) — > 1 as i — > oo 
and also that u( n+1 )(i) and t/ n+1 )(i) both tend to zero as t tends to infinity for 
any n > 1. This shows that = n( n+1 )(t) + uO +1 )(i) -> as t -» oo for 

any n > 1. 

One could hope to deduce further properties of £ by using these integral 
representations. We have not succeeded in doing this, since much cancellation 
between u and v takes place and we do not know A and a explicitly. Further- 
more, it it not even true that + 1) is a Bernstein function even though it is 
increasing and concave. 

We have found an approach using complex analysis and we shall make ex- 
tensive use of the theory of the so-called Pick functions, see (7j. A holomorphic 
function p defined in the upper half plane is a Pick function if it maps the upper 
half plane into the closed upper half plane, or put in another way, if Qp is a 
non-negative harmonic function. A Pick function has an integral representation, 

p{z) = az + b + fc( J ^ z -^Mt), (13) 
where a > 0, b is real and [i is a positive measure on the real line. Furthermore, 

a= i im MM, b = Kp{i) 

y^oo y 

and 

$Sp(x + iy) 

\x = iim ax 

y^o+ 7T 

in the vague topology. Note that a Pick function can be extended to a holo- 
morphic function in C \ R by (|13|) . If A C M is closed we recall that p has a 
holomorphic extension to C \ A if and only if the support of \i is contained in 
A. 

The proofs of the results to follow are based on an investigation of the 
holomorphic function 

f(z) = e~ z + z. (14) 

Concerning the values of the function £ and its derivatives at z = 1 we have 
the following proposition. 

Proposition 5.3 The function ^(z) is holomorphic in a neighbourhood of z = 1 
and we have £(1) = 2/3 and £'(1) = 8/135. 

Proof of Proposition \5.<A The function f(z) — 1 has a zero of multiplicity 2 at 
z = 0. Hence there exists a holomorphic function h in a neighbourhood of 
such that f(z) — 1 = h(z) 2 there. Since h'(0) 2 = 1/2 ^ 0, h is one-to-one near 
z = 0. We choose h f (0) = 1/V2 and with this choice h is uniquely determined. 

We have thus a radius r > such that for any w with \w\ < r, there are 
exactly two solutions to the equation f(z) = w, namely z = h~ 1 (±y/w — 1). In 
particular, for t > 1 and close to 1 we have 

u(t) = hr l (-Vt=l) and v(t) = h^i^/t^T), 
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where u and v are the functions appearing in (JBJ). We denote by X^^Li a nW n 
the Taylor series of h~ 1 (w). Then 

oo 

u(t) + v(t) = ^22a 2k {t- l) fc , 

k=l 

and, since £(t) = u'(t) +v'(t), we find £(1) = 2a 2 and £'(1) = 4a4. The numbers 
a 2 and 04 can be found as follows. If h(z) = Y^=i b n z n then 

\n=l / 

and this gives after some manipulation bi = l/y/2, b 2 = -y/2/12, b 3 = y/2/72 
and 64 = — v2/540. The numbers 02 and 04 can now be identified if we differen- 
tiate four times the relation (h~ 1 )(h(z)) = z and put z = (so that expressions 
involving e.g. (/i _1 )( 4 )(0) appear). We find, again after some manipulation, 
(^~ 1 )W(0) = 2/3 and (/r 1 )^) = 16/45, so that = 2a 2 = 2/3 and 
£'(1) = 4a 4 = 4 • (16/45) • (1/24) = 8/135. □ 

Remark 5.4 One can in principle find any derivative of £ at 1 in this way. 
Computations show that £"(1) = —16/2835. The higher order derivatives of £ 
at 1 can be found in the appendix. 

The monotonicity properties of £ follows from the next proposition. 

Proposition 5.5 The function £ increases on the real line and it is concave 
on [1, 00). 

The key to the proof of this proposition is an integral representation of £ that 
we give in Proposition 15.121 

Lemma 5.6 The function f(z) = e~ z + z is a conformal mapping of the strip 

S 2lT = {0<%z < 2vr} 

onto the domain 

T = C \ = 2iri, dtw > 1} U {Qw = 0, $lw > 1}) . 

In Figure 21 we have indicated the image of some horizontal lines in S 2lT under 
the function /. 

Proof. We have f(x + iy) = a + it if and only if 

e~ x cos y + x = a 
— e~ x siny + y = r. 

For y = 7r we get r = tt and a = x — e~ x so / maps the horizontal line R + in 
onto itself and is one-to-one there. 

Note that for < y < tt we have r = — e~ x siny+y < y so /(<S) C {Qw < tt}. 
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Figure 4: The image of horizontal lines under / 



For r < 7r we get e x = (sin y)/(y — r) and since r < y we find x 
log (sin y/(y — r)) and finally we shall look for a solution y to the equation 



For r < 7r we put 



sin y (y — t) cos y 

log + — - = a. 

y — t sin y 



F r (y) = log + : . (15) 

y — t sin y 



We claim the following: 

1. For < r < ir, F T (y) decreases for y G (r, 7r) from oo to — oo. 

2. For r = 0, ^(y) decreases for y G (r, 7r) from 1 to — oo. 

3. For r < 0, ^(y) decreases for y G (0,7r) from oo to — oo. 

It is easy to see that lim J/ ^ 7r _ F T (y) = — oo, lim^^^ F T (y) = oo for r > and 
lim J/ _ > o+ ^b(y) = 1- We also have lim ?/ _ > o + -^V(y) = oo for r < since the first 
term in F T {y) has a logarithmic singularity. To verify 1., 2. and 3. we need to 
show that F^iy) < 0. We find, after some computation, 

F AV) = 2 cosy : . 

sm y \ y — r sm y / 

If we put k = (sin y)/(y — r) then k > and hence k+1/k > 2. Since cos y < 1 
we see that indeed -F^(y) < 0. 
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From 1. and 3. it follows that for given a G R and r < it, r ^ there is 
a unique solution y to the equation F T (y) = a and therefore there is a unique 
solution to f(x + iy) = a + ir. If r = there is by 2. a unique solution y to the 
equation F T (y) = a when a < 1 and none when a > 1. 

For 7r < y < 2-7T we have siny < so r > y and therefore 

/({vr < 92 < 2vr}) C {vr < 

For t > TT we put 

„ , s , -siny (r-y)cosy 

F T (y) = log + : . 

t — y — sin y 

It follows that F T (y) increases with y E (7r, 27r), and the conformality follows in 
the same way as for the case r < it. □ 

Lemma 5.7 The function w \— * / _1 (1 — w) is a Pick function with the repre- 
sentation 

/-(i - ») = nr'd - o + f - ^) 

w;/iere 

7 ? (t) = i9/- 1 (l-t) 

7T 

is increasing on (0,oo) /rom to 1. 

Proof. First of all, 1 — w belongs to the lower half plane when w belongs to the 
upper half plane. Since / _1 maps all of T, and hence in particular the lower 
half plane, into the strip < < 2tt, / _1 (1 — w) is certainly a Pick function. 
We next derive its integral representation. Since its imaginary part is bounded 
the number a in the representation (|13j) must be zero. 

It remains to identify the measure Since the function / _1 (1 — t — iy) is 
continuous on e.g. R x [0, 1] and is real for y = and t < we find 

Sf-Hl-t-iy) SSf- l (l-t) 
rj(t) = lim — i ^ = — i '- 

for t > and rj(t) = for t < 0. 

By definition the function Y(t) = 7rrj(t) satisfies the equation 

(with the notation of Lemma I5.6() . We find from this 

Y'(t)F^(Y(t)) = -1, 

and since Fq is negative (see again Lemma l5.6j) . Y(t) and hence rj(t) must be 
increasing. Since Y(t) tends to tt as t tends to oo (this is because Fo(Y(t)) — > 
— oo), 7j(t) — ► 1 as i — > cxd. □ 
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Lemma 5.8 The function 

5 H = /- 1 H+7 r V) 

has a holomorphic extension to C \ {Qw = ±2iri, $lw > 1}. 

For t > 1 we have g(t) = u{t) + v(t) and g'(t) = £(t). The function g' is 
thus a holomorphic extension of £. 

All points on the lines {Qw = ±2ni, $lw > 1} are singular points for g. 

Proof. Since / _1 is a conformal mapping of the region 

C \ ({Qw = 2m, $tw > 1} U {Qw = 0,$tw> 1}) , 

the function g is holomorphic in 

C \ ({Qw = ±2ni, > 1} U {9u> = 0, > 1}) . 

However, g has a continuous extension to {Qui = 0, "Stw > 1} with boundary 
values u + v, so from Moreras theorem we conclude that g is holomorphic 
across this half-line. Concerning the boundary values on the half-lines {9w = 
±2iri, > 1} we have (a > 1) 

lim g(a + ir) = lim f~ 1 (a + iT) + f~ 1 (a — 2-iri) 

T — >27T— T — >2lT— 



= v(a) + 2ni + f~ l (a - 2iri) 

and similarly 



lim g(a + it) = u(o~) + 2iri + / 1 (<r — 2iri). 

Therefore g is not continuous across any segment of the half-line 

{Qw = 2m,dtw > 1} 



and so all these points are singular points. Since g{w) = g(w) the same conclu- 
sion holds for the points on the other half-line. □ 



Proposition 5.9 The Taylor series for £(z) centered at 1 has radius of con- 
vergence equal to 2ir and the asymptotic series 



e<*>(i) 

diverges for any z in C 



Proof. We have £ = g r , where g is the function in Lemma 15.81 If the radius of 
convergence of £ at 1 were larger than 2ir then the primitive g would also have 
a holomorhic extension to this larger disk, and this contradicts the fact that g 
has singular points in that disk. 
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Concerning the divergence of the asymptotic series we use that the radius 
of convergence of the Taylor series is finite. It means that 

limsup ^ , 71 > 0, 

fc^oo V / 

and hence that for some e > 0, |^ fc )(l)| 1 / fc /(A;!) 1 / fc > e for infinitely many k. 
From Stirlings formula we have (kl) 1 ^ ~ k/e and therefore 

limsup | 1/k = oo. 

k—*oo 

This shows on the other hand that the asymptotic series diverges for any com- 
plex number z. □ 
In the following log denotes the principal logarithm defined in the cut plane 
C\ (-oo,0j. 

Lemma 5.10 The function 



V(w) = g(logw) = f 1 (logw) + f~ l (logw) 
is a Pick function. It has the representation 

where the function h is given as 

9f/- 1 (logt-»7r) 



h(t) = 1 



7T 



Furthermore, h increases on (0, oo) from to 1. 
For t E [e, oo) we have 

^(r) = ■u(logr) + u(logT). 

Proof. We consider the holomorphic function g from Lemma 15.81 in the strip 
S = {0 < Ssw < it}. Hence the function V(w) = Qg(w) is a harmonic function 
in S. Furthermore, since f~ l maps all of {Qw < it} \ {x \ x > 1} into S, V is 
also bounded there, with the apriori bound — it < V{w) < it. We claim that 
indeed < V(w) < tt for all w G S. 

We consider the boundary values of V. The horizontal line {t + wr, 1 1 G M} 
is mapped by / _1 to itself, whereas {t — in, \ t £ R} is mapped by / _1 to some 
curve inside the strip S. Therefore — itt) has imaginary part greater than 

— 7r and hence the boundary values V(t + in), t G R, are all non negative. Since 
g(w) = g(w), the boundary values V(t),t £ R, are all zero. We conclude that V 
has non-negative boundary values. Since it is bounded, the maximum principle 
in an unbounded region (see e.g. ^Uj) yields that V(w) > for all w £ S. 

Since log maps the upper half plane onto the strip S, the function ty(w) = 
<?(log w) maps the upper half plane into itself, and is hence a Pick function. 
Since is bounded, ^ has an integral representation of the form 

vpH = m® + [°° - -^—] d»(t), 



x yt — W t 2 + 1 
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for some positive measure /x. Here ^^(i) = for t > since V is zero on the 
real line. We thus find the measure \x to be supported on the negative half-line 
with density 

1 — '- L t<0. 

IT 

After making a change of variable (t t— > —t) in the integral we conclude that 

•H-wm-jT^-^) »(*)*. He) 

where /i is the function in the statement of the lemma. 

By definition of h, 7r(l — = 3 ! / _1 (logt — wr), so h is increasing if the 
solution Y(t) to the equation F- w (Y(t)) = logt is decreasing. This is indeed 
the case, since 

F'_^Y{t))Y'(t) = - t > 0, 
and < 0. □ 

Remark 5.11 Since it is easily verified that V(w) = Qg(w) is positive for 
7r < Qui < 2ir it follows that g(2\ogw) is also a Pick function. 

Proposition 5.12 The function ^(logio) is holomorphic in the cut plane C \ 
(—oo,0] with the representation 

f°° t 

e(io gw ) = i- / ——h'(t)dt, 

Jo t + w 

where h(t) is given in Lemma \5.1(A 
If we let w — ► + we see that 

/•oo 

f (log w)->l- / h'(t) dt = l- lim h(t) + h(0) = 0. 
Jo <^°° 

Hence £(t) — > as t — ► — oo. 

Proof of Provosition 15. 1 6 A As we have noticed above, ^(w) = u(logw) + 
v(logw), for w > e , where ^ is the function in Lemma I5.1UI Hence ^(logw) = 
w^f'(w) and so we get from 

f°° whit) , 

If we perform integration by parts on the right-hand side of this relation then 
we get 



or 

coo 



£(logw)= / -—h'(t)dt, 
Jo t + w 

r°° t 

£(lo gU /) = l- / h'{t)dt. 
Jo t + w 



□ 
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Proof of Provosition \5.5\ We get by differentiation of the formula in Propo- 
sition 15.121 that 

m»i = r « mdt , 

w J (t + w) 2 

and since we know that h'(t) > for all t S R, we see that > for all 
t G R. 

It is more technical to establish that £(t) is concave for i > 1. We differen- 
tiate both sides of the relation above once more and we get in this way 

eOogw) - e'Clogw) = ~2W 2 / -—^3 h'(t) dt, 

Jo (t + w ) 

so that, again using the integral representation of £' (log uu), 

C (logw) = / 7 ^ + 7 ttt t/l '(t) dt 

wtn (t) dt 

(* + ™) 3 

f°°s — l 

= w -. r-r sh'(sw) ds. 

Jo ( s + 1) 

This last integral we split into two, one for s < 1 and another for s > 1. We 
thus get, after making the substitution s — ► 1/s in the latter one, 

£"(logw) Z" 1 a-1 ,„ SJ Z" 1 1/s-l ds 
' sh'(sw)ds+ / — yr-h [ — ) 



w 7o (s + 1) 3 ' 7 (l/s + l) 3 s \s/ s 2 



(1 + s) 3 

From this relation we see that £ is concave on [log wq, 00) if 

~ti (—) - sh'(sw) < for s G (0, 1) and to > w . (17) 
S V s / 

Now, as noted in the proof of Lemma |5.1U1 h(s) = 1 — Yl^log s)/tt, where 
Y^{t) £ (0, 7r) is the solution to the equation F_ 7r (y_ 7r (t)) = t, see (fT5|) . There- 
fore /i'(s) = — Yl 7r (logs)/(7rs) and hence 

-h' (—\ - Y -^ logw ~ log,s ) 

S V S / 7TO 

and 

,// x Y^logiu + logs) 
sn (sw) = . 

We see that lfT7)l holds provided 

(log w — log s) > (log u; + log s) . 

Since logs runs through (—00, 0) when s £ (0, 1), this condition is the same as 

r^iogw + t) > YL^iogw-t) 
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for all t > and all w > wq. We aim at proving this for wq = e, so what we 
really should prove is the following: 



f> + t)>r>-t) 



(18) 



for all t > and all a > 1. We verify this inequality by using Lemma |5.7I 
According to that lemma, with w = 1 — a — t + in, we have 



y_ ff (o + t) 



3/ -1 (a + *-***)) 
3f/- x (l- (l-a-i + i7r)) 



7T 



o (s + t + a-l) 2 + vr< 



77(s) ds 



so that 



o 



7T 



^( s + t + a _l)2 + 7r 2 ) ^( S ) ds 







7T 



r»7(s) 



(s + f + a - l) 2 + vr 2 

7T 

o (s + t + a-l) 2 + vr 2 

7T 







/o (s + t + a-l) 2 + vr- 
If we replace t by — t in this formula we get 

7T 



Fl w (a-t) = - / 
J o 



(s - i + a - l) 2 + vr 



t/(s) ds 
»/(s) (is. 



r}'(s)ds, 



so that 

r> + t)-7!>-f) 

/>oo 

4irt 



s + a — 1 



((s + t + a-l) 2 + vr 2 )((s - t + a - l) 2 + ir 2 ) 



rj'(s) ds. 



Since rj'(s) > we see from this expression that Y!_ v (a + t) — Y!_ n (a — t) has 
the same sign as t for a > 1 and hence is positive for f > 0. □ 



6 Relation to work of Ramanujan 

Defining 9{n) for natural numbers n by 



n— 1 i, 
rr 



2 ^ fc! 

fc=0 



??■! 



(19) 



Ramanujan claimed that | < 9(n) < |. This was later proved inde- 
pendently by Szego |12j and Watson ^3]. Further details about Ramanujan's 
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problem can be found in See also [S]. Choi noticed the relation between 
this problem and the median in the gamma distribution, see [5]. By n partial 
integrations we get the formula 



^ ^ Sir 

JU k=0 



hence by (TlT)|) 



and finally 



1 n n 

^e- t dt = --—(l-6(n))e- n , 
n\ n 2 n\ 



/ P \n rm(n+l) 

l-6(n) = \~J J t n e- l dt. (20) 

Following Watson we can write 



n, \ -> e n , n! ^-^ n k 

Bin) = 1 + - — n\ > — 

v ' 2n n n n ^ k\ 

k=0 

„n poo 



1 + / t n e- 1 dt / (n + t) n e- 1 dt 

poo /"OO 

/ t n e n ^dt- / (l + t/n^e^dt. 
Jo Jo 



n 

1+ 2 



We now make the substitution s = 1 + t/n in the last integral and split the 
first integral in two by integrating over (0, 1) and (1, oo). This gives 



/oo 
(ie 1 -*)"^]. (21) 

Watson notices that the right-hand side of this equation makes perfect sence 
when n is replaced by a positive real variable x and the corresponding function 
he denotes y{x). We keep the notation 9{x) and make the substitution t = e~ u 
in the integrals to get 

0(x) = 1 + ^e X (J™ e -*(«+e-«)-« du - J° e -x(u+e-)-u ^ _ 

Introducing the functions u, u from Section we get 

B(x) = 1 + |e x e- xi (e^V^) + e- u(t) «'(f)) dt 

and using 

e-*Wt/(t) + e- u ^u'(t) = £(t) - 2, 
we get after some calculation 



x 



[•OO -1 -1 /'OO 

9{x) = |jf e-^C* + 1) ^ = 3 + 2 / e ~^'(* + !) ! 22 » 
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From (UUJ) we get 

" e(fc)(l) 
fc=0 

We further see that the function 4>(t) from |13l p. 300] is given as 

0(t) = ^(t+l). 

In Watson discusses the behaviour of — 1 _1 log</>(i) for t > 0. He reaches a 
conclusion using a combination of rigorous analysis and numerical tabulation. 
The present analysis does not contribute in making his conclusion rigorous. 

7 Appendix: Higher order expansions and Maple 
code 

In this section we give the derivatives of £ at 1 of order up to 10 and higher 
order asymptotic expansions of (p(x) and m{x). These results we have found 
using the "Maple 9" system. We have also included the Maple code, with a 
short description. 

Computation of the derivatives of £ 

With the notation of Proposition 15.31 we have £^(1) = 2a 2 (fc+i)(& + !)!• Here 
the numbers a k are defined by 

CO 

k=l 

where h" 1 is the inverse to 

oo 

h{z) =J2hz k , 

k=l 

which satisfies (h(z)) 2 = e~ z + z — 1 and h'(0) > 0. 

The procedure h_polynomial (n) returns the list of numbers [bi, . . . , b n ]. These 
numbers are computed by equating the coefficients in the relation 

/ oo \ 2 oo ( k 

\fc=l / k=2 

We get 

l=i 

for k > 2, from which we determine the number b k +i from bi,...,bk- The 
numbers b±, . . . ,b n are just the coefficients in the Taylor expansion of h(z) at 
z = up to order n. We could therefore also have used the command t ay lor 
but we prefer the more direct and efficient method outlined in the procedure 
below. 
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h_polynomial : =proc (n) 
local b, j : 
b[l] :=l/sqrt(2) : 
for j from 1 to (n-1) do 

b[j + l]: = ((-l)-j/(j+2)! 

-sum(b[l+l]*b[j-l+l] , 1=1. . (j-l)))/2/b[l] : 
end do : j : = ' j ' : 
return [seq(b[j], j=l..(n))]; 
end proc; 



The procedure xi_derivatives (n) computes the numbers [£(1), . . . , <^ n_1 Hl)] 
and it uses the procedure h_polynomial. To find these numbers it is enough 
to determine the numbers a±, . . . , a2 n and this we do by equating coefficients in 
the relation h~ l (h(z)) = z. We put 

oo / oo \ k 

l=k \l=l I 

Hence b\^\ = bi for I > 1 and by Cauchy multiplication, 

l-i 

h+l,l = £ bk,mbl-m, 
m=k 

for I > k + 1. Since h~ 1 {h{z)) = z we must furthermore have 

OO / OO \ CO / I \ 

z = E ak ( E 6fc '' z ' ) = E ( E afe6fc >' ) 
fe=l \l=k J 1=1 \k=l J 

Therefore a\b\^ = 1 and 

l 

fe=i 

for Z > 2. Hence a\ = i = 1 /6i - The point is now that a; can be computed 
when we know a\, . . . , a;_i as well as . . . , bij. This is the main idea behind 
the procedure below. 

xi_derivatives : =proc (N) 

local i,p, j ,k,l,temp,a,h,b: 
h:=h_polynomial(2*N) : 
for i from 1 to 2*N do 

for p from 1 to 2*N do 

b[i,p] :=0: 
end do: p:='p' : 
end do : i : = ' i ' : 

for j from 1 to 2*N do b[l, j] :=h[j] : end do: j:='j': 
a[l] :=sqrt(2) : 
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for k from 1 to (2*N-1) do 

for 1 from (k+1) to 2*N do 
temp : =0 : 

for j from 1 to (1-1) do 

temp : =temp+b [k , j ] *h [1- j ] : 
end do : j : =' j ' : 
b[k+l,l] :=temp: 
end do : 1 : = ' 1 ' : 

a [k+1] :=-sum(a[m] *b[m,k+l] ,m=l. .k)/b[k+l,k+l] : 
end do: k:='k' : 

return [seq(2*a[2*j] *(j ! ) , j=l . .N)] ; 
end proc; 

We could also have used the command taylor to compute the coefficients in 
the expansion of h^ 1 (h(z)), but it runs very slowly, so a more efficient program 
is needed. We have found the derivatives of £ to be: 

2 
3 

2 3 
3 3 -5 
2 4 

3 4 • 5 • 7 

2 5 

3 5 • 5 • 7 

2 6 • 281 
3 8 -5 2 -7 - 11 

2 7 • 23 • 227 
3 9 • 5 2 ■ 7 - 11 • 13 

2 8 • 53 • 103 
~3 10 • 5 2 • 7 - 11 • 13 

2 9 • 373 • 4439 • 557 
~ 3 i2 . 54 . 72 . n . 13 . 17 

2 10 • 2650986803 

313 . 54 . 72 . n . 13 . 17 . 19 

2 11 • 6171801683 

314 . 54 . 72 . n . 13 . !7 . ig 

2 12 - 1117-3835213201 
~3 16 -5 2 -7 2 -ll-13-17-19-23' 

The asymptotic expansion of tp 

The asymptotic expansion of (p(x) can be found from relation (|11[). that is 

«*) =eCftt-e t ^ L r°"(- + «- -!)**■+ 

k=0 k=l ' J0 



^1) = 

e'(i) = 

e (2) (i) = 
= 

e (4) (i) = 

e (5) (i) = 

e (6) (i) = 

e (7) (i) = 

e (8) (i) = 

e (9) (i) = 

e (io) (1) = 
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We use two procedures sum_of _int(n,PHI) and xi_sum(n). The procedure 
sum_of _int (n,PHI) takes an expression of the form 



PHI = E 

and computes the sum of integrals 



m 

Ck 



n ( -l\k„k /-PHI 

/ (u + e--l)*du 



M 

k=l 







up to o(x This is done in the same way as in Lemma 14.41 We put 

oo 

( e -« + «-l)*= ^t w u', 
and observe that by Cauchy multiplication 

I _] W 

7i 



_ — (-1) 



i=2 

for ?7i > 2(fc + 1). Since t\ m = (— l) m /m! we can compute the numbers tp. , from 
the numbers ifc-i,- f° r k > 2. The numbers tfc . are then used to approximate 
the integrand (u + e~ u — l) k by its Taylor polynomial centered at u = of order 
n + k. 

sum_of _int : =proc (n,PHI) 

local 1 ,t ,T, s ,k,m, j , temp, tempi : 
for 1 from 2 to (n+1) do 

t[l,l] :=(-l)"l/l! : 
end do : 1 : = ' 1 ' : 

T:=sum(t[l,l]*PHI"(l+l)/(l+l) , 1=2. . (n+1)) : 
s : =-x*T : 

for k from 2 to n do 
T:=0: 

for m from 2*k to (n+k) do 
t [k,m] := 

sum( (-1) "1/1 ! *t [k-1 ,m-l] , 1=2 . . (m-2*k+2) ) : 
T:=T+t [k,m]*PHI-(m+l)/(m+l) : 
end do: m:='m' : 
s:=s+(-l) -k*x"k/k! *T: 
end do: k:='k' : 
temp :=subs(x=l/x, expand (s)) : 
tempi : =0 : 

for j from to (n+1) do 

tempi : =templ+coef f (temp , x , j ) *x~ (- j ) : 
end do : j : = ' j ' : 
return tempi; 
end proc; 
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We notice that the variable s in the procedure sum_of _int(n,PHI) contains 
many terms of order x~ l with / > k. These terms are not needed and they are 
therefore deleted. That happens in the for loop containing tempi. (It improves 
program efficiency) 

The second procedure xi_sum(n+l) simply computes the asymptotic expansion 



end proc; 

To find the expansion of if (x) we combine these two procedures in the procedure 
asymp_phi (n) . This procedure computes the asymptotic expansion of ip(x) up 
to order o(x~ n ~ 1 ). 

asymp_phi : =proc (n) 

local phi,k: 
phi:=l/3/x: 




k=0 



by calling the procedure xi_derivatives (n+1) . 



xi_sum : =proc (n) 



local temp: 

temp:=xi_derivatives(n) : 

return sum (temp [k] /2/x~k, k=l..n); 



for k from 1 to n do 

phi :=xi_sum(k+l)-sum_of_int(k,phi) : 



end do: 
return phi ; 



end proc; 



We have 



10 




k=i 
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where 

1 
3 

29 

3 4 • 5 • 2 
37 

~3 6 - 5-7 
3877 

~3 9 -5 2 -2 2 
8957413 

3 13 -5 3 -7 - 11 
401 • 8842279 

2 • 3 15 • 5 2 • 7 2 • 11 • 13 
356146891 • 2039 
3I8 . 5 4 . 72 . 11 . 13 
216607304027 • 3077479 

~~ 2 3 • 3 21 • 5 6 • 7 3 • 11 • 13 • 17 

31 • 743 • 4569027042343 
323 . 53 . 73 . 11 . 13 • 17 • 19 
71 • 282699240672481 • 1949 • 5113 
2 . 327 . 57 . 73 . 1J2 . 13 . 17 . 19 



The asymptotic expansion of m 

Since m(x) = xe~' fi ^ it is easy to find the asymptotic expansion of m(x) 
when we have the asymptotic expansion of ^p(x). It is done in the procedure 
asymp_m(n), that gives the expansion of m(x) up to order x 1_n . 

asymp_m : =proc (n) 

local temp: 

temp:=subs(x=l/x, asymp_phi (n) ) : 

temp :=convert(taylor(exp (-temp) , x=0, n+2) , polynom) : 
temp :=subs(x=l/x, temp) : 
return expand (x*temp) ; 
end proc; 

We have 

1 9 

m (x) = x - - + > — + o{x ), 

k=l 



Cl = 

C2 = 

C3 = 

c 4 = 

C5 = 

C6 = 

C7 = 

C8 = 

C 9 = 

ClO = 
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where 



mi = 



3 4 -5 
2 3 • 23 



m 2 = 



3 6 • 5 • 7 
2 3 • 281 
3 9 • 5 2 • 7 



m 3 = 



1714 = — 



2 3 • 17 • 139753 
3 13 -5 3 -7 - 11 
2 3 • 708494947 



17l 5 = - 



3 15 • 5 3 • 7 2 • 11 • 13 
2 3 • 140814348739 



m 6 = 



318 . 54 . 72 . n . 13 

2 3 • 7663181003289047 



777,7 = 



321 . 5 6 • 7 3 • 11 • 13 • 17 
2 3 • 653 • 1359581 • 759929 • 3307 



17l 8 = - 



323.56.73.11.13.17.19 

2 3 • 29 • 1376560394479059407 



17lg = - 



3 27 . 5 7 . 73 . 112 . 17 
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